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We investigate the constraints that can be placed on the cosmic string tension by using the current 
Pulsar Timing Array (PTA) limits on the gravitational wave background. We have developed a 
code to compute the spectrum of gravitational waves (GWs) based on the widely accepted one-scale 
model. In its simplest form the one-scale model for cosmic strings allows one to vary: (i) the string 
tension, Gfi/c 2 ; (ii) the size of cosmic string loops relative to the horizon at birth, a; (iii) the 
spectral index of the emission spectrum, q; (iv) the cut-off in the emission spectrum, n„; and (v) the 
intercommutation probability, p. The amplitude and slope of the spectrum in the nHz frequency 
range is very sensitive to these unknown parameters. We have also investigated the impact of more 
complicated scenarios with multiple initial loop sizes a, in particular the 2-a models proposed in 
the literature and a log-normal distribution for a. We have computed the constraint on G/i/c 2 due 
to the limit on a stochastic background of GWs imposed by the European Pulsar Timing Array 
(EPTA). Taking into account all the possible uncertainties in the parameters we find a conservative 
upper limit of Gfi/c 2 < 5.3 x 10 -7 which typically occurs when the loop production scale is close 
to the gravitational backreaction scale, a « FG/i/c 2 . Stronger limits are possible for specific values 
of the parameters which typically correspond to the extremal cases a <C TGfi/c 2 and a S> TG/i/c 2 . 
This limit is less stringent than the previously published limits which are based on cusp emission, an 
approach which does not necessarily model all the possible uncertainties. We discuss the prospects 
for lowering this limit by two orders of magnitude, or even a detection of the GW background, in 
the very near future in the context of the Large European Array for Pulsars (LEAP) and the Square 
Kilometre Array (SKA). 



I. INTRODUCTION 

Cosmic strings are one-dimensional topological defects 
[l| which may have formed in the early Universe during 
various phase transitions expected in Grand Unified The- 
ories (GUTs). Kibble 2] first proposed and investigated 
the production of cosmological scale topological defects 
in the framework of spontaneous symmetry breaking in 
gauge theories. Subsequently, they attracted the interest 
of many cosmologists, because cosmic strings born dur- 
ing a GUT phase transition are a possible source for the 
density fluctuations which eventually led to galaxy for- 
mation The initial enthusiasm diminished some- 
what after the analysis of the Cosmic Background Ex- 
plorer (COBE) satellite data. It was realized that the am- 
plitude of the measured Cosmic Microwave Background 
(CMB) anisotropies on large-scales and the amplitude of 
density fluctuations measured in the galaxy distribution 
on smaller scales could not be reconciled in cosmic string 
models. This ruled them out as the primary source for 
the large-scale structure of the Universe (see, for exam- 
ple, [71,181). These results were confirmed by other ex- 
periments @, but nonetheless, cosmic strings may still 
contribute to the anisotropy seen in the CMB tempera- 
ture, but with less than 10% contribution [10h12| . 

While appearing to be distinctly non-minimal in terms 
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of structure formation - two mechanisms giving rise to 
similar amplitude fluctuations - such ideas are natural 
within a number of well-motivated inflationary models. 
It has been shown that cosmic strings are present in most 
modern inflationary scenarios, be they standard field the- 
ory strings, such as in supersymmetric hybrid inflation 
[l3l [HI, or cosmic superstring models, as in the case of 
brane inflation [15l - ll7j ]. The term cosmic superstrings 
is used to describe cosmic strings which are also funda- 
mental strings but with their tension, or mass per unit 
length, n is reduced from ~ where M p \ is the Planck 
mass, to a lower value, capable of evading constraints on 
the string tension, by a warp factor. See, for example, 

MM- 

Cosmic strings have a wide range of astrophysical sig- 
natures including: ultra high energy cosmic rays [2(| Hl| , 
gamma ray bursts [221 . [23j , radio bursts and synchrotron 
radiation jSS . [25|, Aharonov-Bohm radiation [2(| l27l |. 
gravitational lensing (strong/micro) (28l - [30l |. CMB im- 
prints (non-gaussianity, small/large-scale anisotropies, 
B-mode polarization) [3 Iraq , effects on matter power 
spectra in 21-cm surveys [39|, |4jj. So far no detection 
has been possible, but an interesting opportunity lies in 
their imprint on the stochastic gravitational wave back- 
ground (SGWB) Hll-S|. The SGWB created by a cos- 
mic string network has a very broad spectrum with fre- 
quencies ranging from below the nHz scale to beyond 
the GHz scale, making them a potential source for ev- 
ery present or future GW detection experiment. Such 
broad GW spectra are expected from only some primor- 
dial sources, such as inflation [4!j[5(|, global phase tran- 
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sitions [51| and self-ordering scalar fields [52f . 

A SGWB created by cosmic strings is a possible pri- 
mordial source for detection by Pulsar Timing Arrays 
(PTAs) operating in the nHz frequency range. A PTA 
|53l |54| consists of an ensemble of millisecond pulsars 
(MSPs) that are observed periodically over an extended 
period of time, usually a number of years, which sets the 
frequency probed with maximum sensitivity by the ar- 
ray. The existence of a SGWB changes the Earth-MSP 
distance and therefore it manifests itself as noise in the 
time of arrival of pulses. These timing irregularities will 
have a specific signature which allows us to distinguish 
them from other types of noise 55 1 . The usage of MSPs 
is necessary since this category of neutron stars combines 
a series of characteristics which makes them very stable 
"clocks" , capable of providing high quality timing mea- 
surements. In the case of a SGWB generated by cos- 
mic strings, the quantity we can constrain is their linear 
energy density (or tension in the Nambu-Goto approxi- 
mation), fi, usually expressed through the dimensionless 
quantity G/z/c 2 , where G is Newton's constant and c the 
speed of light. The value of ^ ~ rf 1 in natural units, 
where h = c = 1, is typically related to the energy scale 
of the phase transition, 77, at which they are formed. 

The precise details of the SGWB due to cosmic strings 
are very sensitive to the nature of the string evolution 
and the spectrum of the radiation emitted by the string 
loops that are created due to the intercommutation of 
the long strings. In particular, the spectrum in the band 
probed by PTAs depends on the distribution of loops 
formed and the amount of radiation emitted into high 
frequency modes due to loops formed in the matter era. 
This can be quantified in terms of a spectral index and 
cut-off in the spectrum [lj. Since these details are not 
well understood, this has led to some confusion in the 
literature with a number of inconsistent constraints be- 
ing published, some of them based on the same data [56l — 
l65j | . This is due to the estimates of the amplitude of the 
SWGB being based on different assumptions; some of 
which we believe are too strong for our current level un- 
derstanding of string evolution. In this paper we will use 
the limit on the SGWB imposed by the European Pul- 
sar Timing Array 65] (EPTA) to compute constraints 
on the cosmic string tension as a function of the param- 
eters which describe the details of the string evolution 
and radiation. Taking into account all possible uncer- 
tainties, we define the most reliable, as opposed to the 
tightest, constraint as that which corresponds to the set 
of parameters with the highest upper bound. By param- 
eterizing our ignorance of string evolution we should be 
able to impose a constraint on G[i / c 2 which will never be 
violated. 

The structure of the paper is as follows. In Sec. II 
we give a description of the one-scale model of cosmic 
strings which we used to construct the GW spectra and 
we present a model for the GW emission mechanism from 
cosmic string loops. In Sec. Ill we present the results de- 
scribing the effect of the various cosmic string model pa- 



rameters on the GW power spectrum. The results of mul- 
tiple scale models will also be presented there, in which 
we extend the applicability of the one-scale model to 
more realistic possibilities. Finally, in Sec. IV we present 
robust limits on the cosmic string tension using the re- 
cently published European Pulsar Timing Array limit on 
the SGWB. We conclude with a discussion of the results 
in Sec. V. 



II. MODELLING THE SPECTRUM OF 
GRAVITATIONAL WAVE EMISSION FROM A 
COSMIC STRING NETWORK 

Our calculation of the SGWB expected from a cosmic 
string network is based on the one-scale model [66l - l68| . 
Caldwell and Allen [47| (hereafter, CA92), and DePies 
and Hogan [63| (hereafter, DH07) have previously con- 
sidered this model and our implementation is a combina- 
tion of both their approaches (see also, [6!| for a recent 
investigation). Before going into the details of the one- 
scale model we will give a very brief picture of a cosmic 
string network. 

The basic constituents of a cosmic string network are 
loops and "infinite" (or long) strings; loops so large that 
only a part of them lies within our horizon radius, appear- 
ing as extremely long strings with no ends. These "infi- 
nite" cosmic strings stretch along with the expansion of 
the Universe and oscillate at relativistic speeds. A funda- 
mental dynamical process which impacts the evolution of 
the cosmic string network is intercommutation, whereby 



strings (self-)intersect, exchange partners and form new 
loops [20]. Cosmic string loops have significant tension, 
equal to their linear energy density in the Nambu-Goto 
approximation, so after their formation they start to os- 
cillate relativistically and decay by emitting their energy 
into the "preferred channel" which in the case of local 
strings is thought to be GWs 0> S3 ■ An alternative pic- 
ture is suggested by Abelian-Higgs simulations [7l"l.T72T] 
(see also, [731 ] for a recent discussion), where we have 
the creation of microscopic loops from the string network 
which immediately decay via gauge boson emission. In 
that case, the emission of GWs from the cosmic string 
network is significantly suppressed and the dominant en- 
ergy loss mechanism is field quanta emission, implying 
that the constraints derived here would be invalid. This 
point remains as a caveat to our analysis. 

This mechanism for loop creation and their subsequent 
decay is important since otherwise the cosmic strings 
would dominate the energy density of the Universe rather 
quickly. Instead the loop production allows the network 
to achieve a scaling regime [66] with all the properties 
of the network being related to the horizon radius oc t, 
known as the one-scale model. 

The important features of the one-scale model are: 

i. The evolution of the cosmic string network is con- 
sidered to take place in a homogeneous, flat FRW 
Universe and all the lengthscales of the network are 
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linear multiples of the particle horizon radius. The 
particle horizon radius in such a Universe is 



d H (t) = o(t)t 



dt' 

^0' 



(1) 



where a(£) is the scale factor. The network evolution 
is considered to take place in a co-moving volume of 
size V(t) — a 3 (t)D 3 , where D is an arbitrary length 
scale which will cancel out from our final equations. 



ii. The energy density of the infinite cosmic strings is 



dl(t) 



(2) 



where A quantifies the number of infinite strings 
present in a horizon volume. The value of A was de- 
termined in numerical simulations [74j which suggest 
a value of A ps 52 for the radiation-dominated era 
and A ss 31 for the matter-dominated era (more re- 
cent simulations (75[ give values within 10% of these 
). The A-dominated epoch is believed to be only a 
very small fraction of the cosmic history, so we used 
the same value for A as in the matter-dominated era 
presuming that our results will not be radically dif- 
ferent in the case of the SGWB calculation. 

iii. Every cosmic string loop is born with a size £b(tb) 
which is a constant fraction of the horizon radius at 
the birth time th 



4(*b) = adB_(t h ) 



(3) 



A correction has to be applied in Eq. © because 
cosmic string loops are born with relativistic peculiar 
velocities which are very quickly reduced due to the 
expansion of the Universe, resulting in an energy loss 
from the loops almost immediately after their birth. 
By taking into account this redshifting of their initial 
velocity we can write the initial length (energy) of a 
newborn cosmic string loop as ^b(*b) = fr<%dn(%) , 
where / r incorporates this energy reduction and is 

fr « 0.71 0. 



A. Number density n(£, t) of cosmic string loops 

In the one-scale model the "infinite" string network 
needs to lose sufficient energy so as to maintain the scal- 
ing regime which allows us to calculate n(£, t)di, the num- 
ber density of cosmic string loops with lengths between I 
and I + dl at time t after the creation of the cosmic string 
network. This quantity gives us a full "map" of the num- 
ber and size of cosmic string loops present throughout the 
cosmic history. In order to calculate the number density 
of cosmic string loops we need information about their 
birth rate. Allen and Caldwell (4]],[7Zl calculated in detail 
the formation rate of cosmic string loops from the conser- 
vation of the cosmic string stress-energy tensor. Briefly, 



they did this by first calculating the equation of energy 
conservation of the cosmic string network 



dt 



[a 3 (t)(p tot {t)+p tot (t))] = p tot (t)a 3 (t) 



(4) 



where ptot and ptot are, respectively, the energy density 
and pressure of the whole system (infinite cosmic strings, 
cosmic string loops and GWs emitted). Combining this 
with the equation of state of infinite strings 



Poo{t) = \ Paa {t) [2(v 2 )/c 2 - 1] 



(5) 



where (v 2 ) is the mean squared velocity of infinite cosmic 
strings, they derived the amount of energy lost by the 
network to create new cosmic string loops per unit time. 



dE; 



loop,cr 



dt 



-V(t) 



Poo(t) 



a(t) 
a(t) 



,{t)(\ + (v 2 )/c 2 ) 



(6) 

The mean squared velocity of infinite strings is also deter- 
mined by numerical simulations [lH and is (v 2 ) /c 2 = 0.43 
for the radiation-dominated era and (v 2 )/c 2 — 0.37 for 
the matter-dominated era. The latter value will also be 
used in our calculations for the A-dominated era. The 
most recent evolution simulations [75j suggest similar val- 
ues. In the one-scale model, since we know the size of the 
newborn loops, we can write 



loop,cr 



dt 



fiadn(t)c 



dN, 



loop 



dt 



(7) 



where N\ oop is the total number of loops created since 
the creation of the network within the volume V(t) and 
dN\ oop /dt is the corresponding formation rate. Combin- 
ing Eqs. ©,© we get 



dN] 



loop 



V{t) 



dt 



[iadn(t)c 2 



,(t) + 2^ 00 (t)(l 
a(t) 



2 )/c 2 ) 



(8) 



and by using Eq. © we can bring it to the simpler form 



rfMoop = 2V(t) Poo (t) 
dt [iadn(t)c 2 



d{t)(v 2 



d H (t) a{t)c 2 



(9) 



Knowing the formation rate, we can calculate the number 
of loops born at any instant in cosmic history. Note, that 
this rate was actually calculated from the energy lost by 
the network in order to maintain scaling, which means 
that our results are automatically normalized. This is 
not always the case in other methods. 

Each cosmic string loop decays by GW emission with 
a constant rate 



dE 



loop, em 



dt 



-YGpfc, 



(10) 



where E\ ooPiem is the energy emitted by a cosmic string 
loop in the form of GWs and T is a constant which de- 
scribes the efficiency of the emission mechanism. The 
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value of r depends on the shape of the cosmic string 
loops. Throughout this work, a value of T w 50 will be 
used which has been calculated as the average for repre- 
sentative loops through numerical simulations [Z^L An- 
alytic calculations of T have also been performed [79l l8fjj 
but only for specific cases. Assuming that GW emission 
is the dominant energy loss mechanism of cosmic string 
loops, the length £(t,t\,) of a cosmic string loop at time 
t, if this loop was born at time th, can be written as 



£(t,t h ) = / r ad H (*b) 



TGfi 



(t - t b ) . 



(11) 



This relation for the variation length of a cosmic string 
loop with the cosmic time and its birth time can be in- 
terpreted in an inverse way; a cosmic string loop which 
has a length equal to £ at time t, has to have been born 
at a specific time t\, given by the solution of Eq. (TTTT) . 
and only at that time. 

We will express the n(£,t) function as a discrete two 
dimensional array, with cosmic time and cosmic string 
loop length as the axes. Each element n(£i,tj) of this 
array we will be calculated in the following way. First, we 
calculate the corresponding birth time, tbj , of these loops 
from Eq. (TTTI) . We define Af(£,t)d£ to be the number of 
loops with length between £ and £ + d£ present in our 
volume at time t. Of course, Af(£i, tj) = Af(£b,j, ib,j)- For 
the number of loops born at time tb,j in our simulation 
volume we will have 



Af(£,t)d£\ 



t=t h j,i=i b 



. — dNi oop \ . 
,j fit— i hi 



(12) 



Since d£ = (/raring) — TGfi/cjdt, by substituting in 
Eq. (Tj"2"|) and dividing by our simulation volume V(ij) we 
find that 



n(£ u t 3 ) = 



1 



loop 



/ r ad H (tb,j) + TGfi/c 



dt 



(13) 

This allows us to calculate all the elements of the n(£, t) 
array and, by using multi-variate interpolation, we can 
construct the n(£, t) function. Of course, we do not have 
to calculate the elements for which £\ > f T adti(tj) which 
are equal to zero. The veracity of our results was tested 
against the analytic formulas for the loop number density 
in the radiation and matter dominated eras [76] for a 
range of values of a and TG/j,. The more recent and 
detailed analytic formulas found in [8l[ yield expressions 
which have to be fitted to evolution simulation results 
and do not have the flexibility of those by Bennett and 
Bouchet. 



B. Gravitational wave emission from string loops 

A cosmic string loop oscillates relativistically under its 
tension and emits GWs in a series of harmonics with 
frequencies which depend only on the length of the loop, 
£, and the harmonic mode n. The period of the loop is 



£/2c and the particular frequencies are harmonics of £/2, 
that is 



_ 2nc 



(14) 



Given a particular string trajectory, in the Nambu-Goto 
approximation and ignoring the effects of radiation back- 
reaction, one can compute the total radiated power from 
a loop, P — TGfi 2 c, which is independent of i. The 
value of r depends on the specific trajectory, but as al- 
ready noted it is s» 50 for typical loop trajectories [78j . 
The power emitted into each harmonic mode is given by 



dE, 



gw,loop 

~dt 



= P n G f i 2 c, 



where 



P n = Tn-y J2 



(15) 



(16) 



is a coefficient for each mode which determines the 
amount of radiated energy that is emitted through the 
respective mode and q is the spectral index. The value of 
q can be computed for a specific trajectory: q = 4/3 for 
the Kibble- Turok loops and q = 2 for a square loop with 
kinks [l|. It can be argued that q 4/3 for any string 
loop which has a cusp; something which is expected for 
string trajectories without a kink. 

The effects of radiation backreaction are ignored in this 
calculation. The decay of the loop length can be de- 
scribed by the linear decay with time already discussed in 
the previous section. However, the precise details of the 
spectrum can have significant impact on the amplitude 
and slope of the SGWB in the nHz region that is relevant 
to PTAs and this can also be affected by backreaction. 
It has been shown in full field theory simulations that 
the spectrum of Goldstone boson radiation from global 
strings is significantly softened by the effects of radia- 
tion backreaction [821 ] . In particular, it was shown that 
an initial q = 1 spectrum was dominated by the funda- 
mental mode after a small number of oscillations. The 
equivalent simulations are not possible in the context of 
gravitational radiation, but there are sufficient similar- 
ities between the radiative mechanisms to suggest that 
something similar may also take place in this case as well. 

In order to model this effect we include an extra phe- 
nomenological parameter, n», as was done in (48[: P n = 
for n > n* and the normalization factor is modified for 
n < n* so that 



Pn = Tn-y 



m— 1 



(17) 



and the total power emitted is unchanged. The value 
of n* is unknown, but it must be less than the ratio of 
the loop length to the string core width, 5. For macro- 
scopic strings, £/ 5 is typically very large, but the round- 
ing of the cusps and kinks expected in a realistic network 
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might mean that n* « i?/<5 where i? is the local radius 
of curvature which could be much less, say in the range 
log 10 n* = 3 to 5. In what follows we will allow and q 
to be free parameters which wc will vary. 



C. Stochastic gravitational wave background 

The standard quantity used to quantify the amplitude 
of the SGWB is energy density in GWs per logarithmic 
frequency interval measured relative to the critical den- 
sity, p cr it, which is given by 

1 dp„ m 



^gw(/) 



p cr it a log / 

This can be related to the dimensionless strain of the GW 
by 



h gw (f) = 1.3 x 10- 9 ^ gw (/)^ ( 



InHzN 



(19) 



Since p C rit = 3Hq /8nG depends on the Hubble Constant 
Hq = 100/ikms -1 Mpc -1 , it is conventional to plot the 
dimensionless quantity f2 gw (/)/i 2 . 

The starting point for our calculation is based on the 
formula that gives the spectral density of the emitted 
GWs derived in PJ, 



rf r ad H (t') 

X / 

JO 



td£n(e,t')g(^^fl \ , (20) 

where / is the frequency of the GWs as we observe them 
today, t s» 13.4 Gyr is the present time, tf = i p ic 4 /(G/i) 2 
is the time of formation of the cosmic string network 
[47l ] with i p i the Planck time and g(z) a function which 
describes the spectrum of radiation emitted by a loop 
and is normalized by J Q g(z)dz = YGp?c. We note that 
in [l[ the integral is written in terms of angular frequency 
uj = 2nf and g(z) is normalized by T with the factor of 
G/i 2 c included in the equivalent of Eq. (l20l) . 

We will model a discrete emission spectrum and there- 
fore g(z) will be a sum of ^-functions given by 



g(z) = Gp 2 c J^PjSiz-Anj). 



(21) 



(22) 



If we set z = (a(t )/a(t'))(2Trf£/c) then 

and hence substituting Eq. (|2"2"j) and Eq. (|2ip into 
Eq. ([H we find that 

dp sw 2Gp 2 c 3 r t0 



df P 



dt 1 



A \ 5 n 



Y.i p A a ^ 2 4A- (23) 



°(*o)/ hi W*o) f 



The integral of Eq. (|2"3"|) requires the continuous calcu- 
lation of the argument of n(£, t) during its numerical eval- 
uation. We can decouple these two calculations, some- 
thing which also gives a better physical intuition into 
what Eq. (|2"3")) represents. A similar approach was also 
followed in DH07, but note that some quantities are ex- 
pressed differently in our implementation. As we have 
already mentioned, the sum in Eq. (1231) gives the contri- 
bution to the GW spectrum of each emission mode from 
each cosmic string loop. The quantity a(t')c(a(to)f)~ x 
in the number density is actually the length of the cos- 
mic string loops which emit GWs at time t' which are 
observed at the present day with redshifted frequency /. 
This loop population is the only one we are interested 
in when we have to evaluate Eq. (|23p in a specific fre- 
quency bin. It is reasonable then to express the integral 
in terms of a new function, say n(/, t), which gives the 
number density of loops which at time t emit GWs ob- 
served today with frequency /. In this case, f2 gw (/) can 
be written as 



^gw(/) 



2Gp 2 c 



2„3 



\ ' 



Pcrita 5 (i )/ 2 j-[ jtt 



III 



jPj a 5 (t') nj (f,t')dt' . 



(24) 

Of course, the n?(/, t) function has to be constructed for 
every emission mode, j. 

The rij (/, t) function can be constructed in a straight- 
forward way once we have obtained n{i,t). Using the 
same array we constructed for n(£,t), we change the cos- 
mic string loop length axis with the corresponding fre- 
quency axis by converting all lengths to frequencies using 
Eq. (1141) . Multi-variate interpolation with this new axis 
will give the rij(f- ln .,t) function which gives the number 
density of loops which at time t emit GWs at frequency 
/in.. We will use this rij(f iD ,t) to construct n,j(f,t). We 
can calculate the elements nj(f m ,t n ) of the rij(f,t) data 
array in the following way. First, we calculate the fre- 
quency / mj i n . which the GW had when they were emitted, 
/m.in. = (a(to)/a(t n ))f m . Once we have this information, 
we use the previously constructed n(/ n .,i) function and 
calculate the respective loop number density. In this way 
we populate the n(f, t) data array and we create the cor- 
responding interpolating function. 

A typical GW spectrum from cosmic strings is pre- 
sented in Fig. [TJ It has two distinct features: a low fre- 
quency peak and a flat spectrum at higher frequencies. 
The flat part of the spectrum is created from GWs emit- 
ted by cosmic string loops which were created during the 
radiation era. Roughly, it is flat from > 1 nHz to around 
the GHz frequency range and then it drops rapidly in 
amplitude. The low frequency peak, which dominates 
below ~ 1 nHz, is created from the more recent emission 
of cosmic string loops formed during the matter era. In- 
terestingly the cross-over region between the two regimes 
is exactly that which is probed by the PTAs. As we will 
see, the point at which this cross-over actually takes place 
is sensitive to the parameters describing the string net- 
work and the radiation spectrum. 
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FIG. 1: The GW energy density per logarithmic frequency 
interval n gw (/)/i 2 of a cosmic string network with Gp/c 2 = 
1CT 7 , a = 10~ 3 and n* = 1. The black (solid) line is the full 
spectrum from the network due to loops formed in both radi- 
ation and matter eras, whereas the red (dashed) line is that 
from the radiation-dominated era and the blue (dot-dashed) 
line is from the matter-dominated era. The grey shaded area 
shows the frequency window probed with the highest sensi- 
tivity by PTA experiments with duration between 5 and 10 
years. 



D. Intercommutation probability 

Whenever two field theory cosmic strings collide they 
exchange partners with an intercommutation probability 
p = 1 [70] . This is not necessarily the case for cosmic su- 
perstrings however, which intercommute with a reduced 
intercommutation probability p < 1. This can be at- 
tributed to the extra dimensions in which cosmic super- 
strings are moving, with a successful intercommutation 
requiring their collision in all dimensions and not just in 
the three spatial dimensions visible to us. If p < 1 then 
the scaling density of long strings is increased in order 
to increase the number of intersections per unit time and 
hence allow the network to lose the requisite amount of 
energy necessary to maintain scaling. This will increase 
the number of loops and hence will increase the ampli- 
tude of the SGWB by a uniform scaling. There is, how- 
ever, some controversy as to the exact dependence on p. 
Jones, Stoica and Tye [l9j . argued that the self-similar 
length scale, L, of the cosmic string network should scale 
as L cx pt, which would mean that poo oc L~ 2 oc p~ 2 . 
In that case, even a small decrease in p would lead to a 
dramatic increase in the amplitude of the SGWB. How- 
ever, in such a case the inter-string distance d s , due to 
the higher string density, is smaller than the length scale 
of the network L, whereas in the one-scale model L ~ d s , 
suggesting that this argument needs to be modified. 

Sakcllariadou [83[ has performed simulations of cosmic 
superstring networks in Minkowski spacetime which sug- 
gest that L oc p x l 2 t, implying that poo oc It was 
suggested the discrepancy with the results of Jones et 
al. stems from the small-scale structure of cosmic stings, 



which ensures more intersection points when two strings 
collide, and therefore there are more chances for success- 
ful loop production. 

There are two techniques used to model the dynam- 
ics of strings in the Nambu-Goto approximation: one 
is the Minkowski spacetime approach used in (83|; the 
other is to model the expansion of the Universe. The 
results of such simulations are reported by Avgoustidis 
and Shellard in 0, HI- They find that when p < 0.1 
then poo cx p~ 6 , whereas for 0.1 < p < 1.0 they find 
Poo a They also suggest that small-scale structure 
is responsible for the difference from the poo oc p~ 2 scal- 
ing law and they propose a simple two-scale model which 
describes quite accurately their simulation results. The 
difference in the scaling laws of [83[ and [85| has to do 
with fitting model parameters to results of fundamen- 
tally different simulations, so the exact reasons for this 
discrepancy are not easy to trace. 

In this work we will not make a judgement on the pre- 
cise dependence of the scaling density of infinite strings 
as a function of p except that it can be modeled by a 
power law 



A( P ) 



A(l) 



p" 



(25) 



where k is the model parameter and A(l) = 52 and 
A(l) = 31 in the radiation and matter eras respectively. 
The results of [83| suggest that k = 1, whereas those 
of [H, HI suggest k = 0.6 for p < 0.1 and k = 1 
for 0.1 < p < 1.0. The consequence of this assump- 
tion is that the amplitude of the SWGB will scale as 
O gw (/) cx p~ k independent of /. 



III. 



CHARACTERISTICS OF COSMIC STRING 
INDUCED SPECTRA 

Low frequency cut-off due to newborn large 
loops. 



As we mentioned in Sec. Ill Bt each cosmic string loop 
emits GWs into an ensemble of harmonics defined by 
/„ = Incjl. This means that there is a low frequency 
cut-off on the GWs that a cosmic string network emits, 
defined by the first emission mode of the largest loops 
present. The largest loops are those created at the 
present time to and have length £g — f r adn(to), with 
a corresponding low frequency cut-off /o cx l/ato- The 
redshifted frequencies of the GWs emitted by loops pre- 
viously born will always be higher than /o in both the 
radiation- and matter-dominated eras. For example, in 
the radiation era the frequency of the first emission mode 
of a loop formed at time t\ redshifted to the present is 

/i oc te/ 6 /a(*i /2t o /3 ) > fo, where t cq a 25,000 yrs is the 
time of radiation-matter equality. The same calculation 

in the matter era gives /i oc l/a(tY 3 t^ 3 ), which is also 
greater than Jq. To demonstrate the strength of this in- 
equality, in the matter era, the GWs of the first emission 
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FIG. 2: The GW sensitivity curves for a 10-year (black thick 
line) and a 5-year (black dashed line) PTA experiment, with 
the 10-year experiment achieving slightly better maximum 
sensitivity. The frequencies where these experiments achieve 
maximum sensitivity are 3.2 nHz and 6.3 nHz respectively. 
The red thick line is the GW spectrum of a cosmic string 
network for ai = 5.7 x 10 -10 and the red dashed line is the 
spectrum for a-z = 2.8 x 10 -10 network. While the 10-yr ex- 
periment has a greater overall sensitivity at its minimum fre- 
quency, it has a lower sensitivity at the frequencies to which 
the 5-year experiment is sensitive to (see text for details). 



mode of a loop with a = 0.1 emitted at the time of its 
birth, say t\ = 10 10 Gyrs, will be observed today at a fre- 
quency Ri 2.7xlO~ 17 Hz, whereas the corresponding GWs 
of a loop born at the present time will have a frequency 
ss 1.8 x 10~ 17 Hz. Similarly, the GWs of a similar loop 
born in the radiation era at time, say t\ = 10 4 yrs will 
have at the present time a frequency « 4.3 x 10~ n ,Hz. 

This low frequency cut-off needs to be treated very 
carefully because of its strong dependence on a, which is 
unknown. For small, but totally acceptable values of a 
(such as, a ~ 10~ 12 — 10~ 16 ), the low frequency cut-off 
can be as high as a few microhertz. Although, ground- 
based observatories and LISA are not seriously affected 
by this low frequency limit since they are sensitive to 
frequencies ~ 100 Hz and ~ 10~ 3 Hz respectively, PTAs 
can easily be rendered useless for detecting emission from 
such cosmic string networks. As we will show in the 
following paragraphs, PTAs can adequately probe the 
GW emission of cosmic string networks with a < 10~ 9 . 

A PTA sensitivity curve roughly has the shape of an 
inverted triangle, with its sensitivity peaking at wave- 
lengths similar to the duration of the PTA experiment. 
In Fig. [2] we plot the sensitivity curves for two completely 
different PTA experiments, one with 5 and one with 10 
years duration, which achieve maximum sensitivity at 
frequencies ~ 6.3 nHz and ~ 3.2 nHz respectively. For 
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FIG. 3: Regions of the a — n parameter space which can be 
probed by PTA experiments. The dark gray region includes 
all the cosmic string network configurations which create a 
SGWB probed at maximum sensitivity by a 10-year PTA ex- 
periment. Additionally to this region, the light gray slice 
includes all the extra configurations which can be probed at 
maximum sensitivity by a 5-year PTA experiment. The white 
area includes all those configurations which are probed by the 
reduced sensitivity slope (see, Fig. for both 5- and 10-year 
experiments. The hatched area includes the configurations 
which are inaccessible to PTAs. 



frequencies higher than these, the sensitivity decreases 
with a slope given by Rf, where R is the root- mean- 
square of the residuals in pulsar timing data [86| and / 
the frequency of the GWs. For the sensitivity curves of 
Fig. [2] we assumed that the maximum sensitivity of the 
5-year experiment is at fi gw /i 2 = 1.2 x 10~ 8 and that 
the 10-year experiment has a slightly better sensitivity, 
equal to fi gw /i 2 = 8.9 x 10~ 9 . This assumption actually 
implies that the RMS residuals in the 10-year PTA ex- 
periment are larger than those in the 5-year experiment. 
The RMS residuals of the time of arrival of pulses are 
expected to improve as oc N^ 1 / 2 , where N is the num- 
ber of the time of arrivals (observations) used. There- 
fore, if the 10-year PTA experiment was performed at 
the same telescope/instrumentation with the 5-year ex- 
periment we would expect a much more improved sensi- 
tivity on Q gw h 2 , due to the double amount of data. The 
value of the highest frequency in the sensitivity curve 
depends on the mean time between observations of the 
same MSP in the whole data span. As an example, for 
the EPTA where each MSP is typically observed once 
every two weeks, the maximum freq uency is ~ 830 nHz. 
The PPTA [13 and NANOGrav follow a similar ob- 
serving schedule. 

Now, let us consider the sensitivity of PTAs to emis- 
sion from string loops with size ~ at emitting into the 
nth harmonic. In Fig. [3] we present the regions in the 
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a—n parameter space which can, or cannot, be probed by 
present and future PTA experiments. The dark gray area 
includes all the a — n combinations which give a SGWB 
with a low frequency cut-off lower than the frequency at 
which a 10- year PTA experiment achieves its highest sen- 
sitivity. A shorter, 5-year experiment is sensitive to all 
these cosmic string networks plus the networks included 
in the light gray slice. Although counter-intuitive, by 
increasing the duration of a PTA experiment we reduce 
the number of different cosmic string network configura- 
tions which are observable at maximum sensitivity, even 
though we increase its overall sensitivity by collecting 
more data. 

We can easily demonstrate this in Fig. [3J where we 
present the GW spectra for two cosmic string networks 
with ol\ — 5.7 x 10~ 10 (thick red curve) and oti = 
2.8 x 10 -10 (dashed red curve) respectively. Both net- 
works have G/i/c 2 = 1CP 7 , n.* = 1 and q = 4/3. The 
specific values for a\, «2 were selected so that their low 
frequency cut-offs coincide with the frequencies which a 
5-year and a 10-year PTA experiment achieves maximum 
sensitivity. In this example, the 5-year experiment probes 
part of both spectra, while the 10-year experiment probes 
the o;i network and just misses the a-i network. It is 
clear, that any network which has a low frequency cut- 
off between 3.2 nHz and 6.3 nHz will be probed with the 
reduced sensitivity slope from the 10-year experiment, 
something that doesn't happen with the 5-year experi- 
ment. It is interesting to sec, that in the specific exam- 
ple, the 10-year experiment will not be able to detect the 
emission from the a.2 network at all. Although in real- 
ity a 10-year experiment is expected to have much better 
RMS residuals of a 5-year one, such an event might be 
true when comparing PTA experiments of different ob- 
servatories. The white area of Fig.[3]contains all the a — n 
combinations which can be probed by both experiments 
with reduced sensitivity, and therefore, higher G/i/c 2 val- 
ues are required to make a detection. The hatched area 
corresponds to the a — n combinations inaccessible by 
PTAs due to their high emission frequencies. From this 
plot we can see that networks which produce very small 
loops with a < 10~ 12 cannot be detected by PTAs, re- 
gardless of their tension. 



B. Effects of cosmic string model parameters on 
the gravitational wave spectrum 

Here, we present a detailed analysis of the effects of 
each cosmic string model parameter on the GW spec- 
trum with the objective of building up a picture of how 
the various uncertainties can affect the spectrum. We will 
often make the distinction between whether loops at for- 
mation are smaller, a < TGfi/c 2 , or larger, a > TGfi/c 2 , 
than the gravitational backreaction scale; we will refer 
to these as small and large loops respectively. Since the 
time of death of the loop, id, and its time of birth, %, are 
related by t& ~ [1+ (cue 2 /TGfi)]tb, this distinction corre- 



sponds to the loops either dying within a Hubble time, 
or living for much longer, respectively. The fiducial set 
of parameters is G/t/c 2 = 10~ 7 , a = 10~ 7 , q = 4/3, 
n* = 1, p = 1 and in the subsequent discussion we will 
vary each of them while keeping the others equal to their 
fiducial values. In our computations we used a scale fac- 
tor suitable for a radiation-matter-A Universe. For the 
numerical values of the constants entering in the com- 
putation of the scale factor we used the WMAP 7-year 
results [HI and in particular the WMAP+BAO+i? pa- 
rameter values, h = 0.704 and Q\ = 0.728. We also used 
f2 r /i 2 = 2.47 x 10~ 5 and Sl m = 1 - tt A - f2 r for a flat 
Universe. 



1. Varying G/i/c 2 

One would expect heavier strings, with larger values 
of G/i/c 2 , to lead to larger SGWB amplitudes and in- 
deed this is typically the case. However, when varying 
G/i/c 2 there are other more subtle effects which can have 
an impact on the GW spectrum. We present plots of 
f2 gw /i 2 in Fig. 2] for different values of G/i/c 2 keeping a, 
q, n* and p the same. In the small loops regime (thin red 
lines), lowering the string tension reduces the amplitude 
of the spectrum with the shape remaining the same; one 
finds that f2 gw /i 2 oc G/i/c 2 . However, in the large loops 
regime (thick blue lines), along with the expected de- 
crease in amplitude, there are also changes in the shape 
of the spectrum. The peak frequency starts to shift to- 
wards higher frequencies, initially being oc (G/i/c 2 ) -1 / 4 
and rather quickly it settles to being oc (G/i/c 2 ) _1 . The 
SGWB amplitude also decreases oc (G/i/c 2 ) 1 / 2 , slower 
than in the case where a < TGfx/c 2 , which can clearly 
be seen in the flat part of the spectrum. 

The frequency of the peak of the spectrum can be ap- 
proximated analytically. The key question in making 
such an approximation is which loop population is re- 
sponsible for the emission at the peak frequency. To 
answer this question, we need to define the birth time 
of loops tb(t) as being the time of birth of a loop which 
dies at time t. The birth time of loops which die at the 
present time, assuming that we are in the matter era, can 
be calculated from Eq. (TTTT) setting I = 0, 

*b(*o)=( 1 + ^^) ^o- (26) 

In CA92, they suggested that the peak emission is cre- 
ated by the n = 1 emission of the most recent, "domi- 
nant population" of loops. Since the birth rate of loops 
is continuously decreasing with time (oc t~ 2 ), they as- 
sumed that this "dominant" population was born at a 
time ~ 2tb(to). This leads to a simple approximation of 
the peak frequency given by / p0 ak = 2c 2 /TG^ito. This 
approximation is presented in Fig. 0] with a long dashed 
red line. It is obvious that since this equation is inde- 
pendent of a it will not give a correct description on the 
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FIG. 4: Plots of normalized gravitational wave energy den- 
sity per logarithmic frequency interval, Q gw h 2 , due to cosmic 
string networks with different tensions but the same fiducial 
values of a, n*, q and p. The thick blue lines are for networks 
in the large loop regime and the thin red lines are from net- 
works in small loop regime. The dashed black line signifies 
the network for which a = TG/i/c 2 . The analytic approxima- 
tions of the peak frequency are also shown: the approximation 
found in CA92 (red long dashed curve) and our improved ap- 
proximation (short dashed green curve). 



different behavior in the large and small loops regimes 
and indeed this can be seen in Fig. @] Moreover, even in 
the large loop regime where it seems to be in reasonable 
agreement, the more we decrease the string tension the 
worse the approximation becomes. 

We have managed to construct a better approximate 
formula for /peak, where we do not make any assumption 
about the birth time of the loop population responsible 
for the peak emission. Instead, we created a general, 
approximate formula and we determine when these loops 
were formed by comparing the analytic results with those 
of our computations. 

The peak frequency must originate from the redshifted 
emission in the n — 1 mode of this population, the lowest 
frequency it ever emitted. Using Eq. (I26[) for the birth 
time of loops we introduce the concept of loop genera- 
tions, g. We will refer to loops which die right now, and 
therefore, were born at time t\ — tb(io)j &s generation 
g = 1 loops. The loops of generation g = 2 are those 
which died when the loops of g = 1 were born and have 
a birth time ti = t\,{ti). In the same way, the loops of 
generation g are those which die when the loops of gen- 
eration g — 1 were born. From Eq. (126[) we find the birth 
time t g of generation g loops to be 



tg = I 1 



3/ r ac 2 



to. 



(27) 



The lowest GW frequency (n — I) emitted by loops of 
generation g in the matter era is 
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FIG. 5: f2 gw /i 2 for cosmic string networks with different values 
of a and the fiducial values of G/x/c 2 , n*, q and p. With thick 
blue lines we plot the networks in the regime of large loops and 
with thin red lines the networks in the regime of small loops. 
With dashed line we plot the network with a — FGfi/c 2 which 
signifies the critical point after which we have no amplitude 
decrease. 



and when we redshift it to the present day, its observed 
frequency is 



fg 



a(t g ) 



a (t ) 3/ r at g 3/ r at 



3/r«C 2 



9/3 



(29) 



Eq. (|29p is the general approximation for the peak fre- 
quency, without making any assumptions about which 
generation's loops created it. Using the results of our 
computations, we found out that the best approximation 
to the peak frequency is given by 



Xpcak 



3/rCrfo 



3/ r ac 2 
TGfj, 



10/9 



(30) 



which is plotted with a short dashed green line in Fig. [4] 
This means that the peak region is due to loops of gen- 
eration g ~ 10/3, i.e. of loops born just before the third 
generation loops. We have changed the numerical factor 
in the parenthesis of Eqs. (|29|) . f|30[) from 1 to 2, so to 
achieve a perfect fit. In any case, this is a minor cor- 
rection (less than 3%) which only affects networks with 
TG/i/c 2 > a. 



2. Varying a 

The effects of varying a in the large/small loop re- 
gions are the inverse of those seen when varying G/x/c 2 . 
In Fig. [S] we present the GW spectra for cosmic string 



networks with the fiducial values of Gfi/c 2 



q and p 



for various values of a. 

In the large loop regime (blue thick lines), as a de- 
creases the most prominent feature is a decrease of the 
amplitude of the overall spectrum. This decrease is 
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FIG. 6: figw/i 2 for cosmic string networks with the same fiducial string tension value, but with different n, and q values in the 
cases of large (a = 0.1) and small (a = 10 -9 ) loops. Different colors are used for different values of n», with the specific values 
denoted in each plot. In the q — 2 case, the spectra for n* = 10 3 are plotted with a dashed line instead of a different color to 
better distinguish it from the n» = 10 2 case since the results are almost identical for extended frequency ranges. 



oc a 1 / 2 when a ^> TG/x/c 2 , but the dependence becomes 
weaker, being oc a 1 / 4 when a gets close to the critical 
value, a — TGfi/c 2 . The higher SGWB amplitude for 
large values of a is expected: large a means the loops per- 
sist for longer periods of time, and therefore have more 
time to emit their energy as GWs. Our results in this 
regime agree with those in DH07. 



The situation is very different in the small loop regime 
(thin red lines) . There is no decrease in the overall ampli- 
tude, nor a significant change in the amplitude difference 
between the peak and the flat part of the spectrum once 
a < TG/i/c 2 . Instead of this we see a shifting of the 
spectrum to higher frequencies, something which agrees 
with the results presented in [62|. The overall ampli- 
tude invariance to changes in a is a radiation era effect, 
where the small loops decay in less than a Hubble time. 
The independence of fl gw h 2 from the value of a can be 
clearly seen in the analytic approximations for ft gw h 2 in 
the radiation era (i.e., see equation 4 in 12]) if we as- 
sume a <?C TG/i. The shifting of the peak frequency is 
oc a -1 , consistent with the minimum frequency emitted 
by a loop being / m i n w 2c/{att > ) and Eq. ([50]) . 



3. Varying the emission spectrum parameters q and n* 



The spectrum of gravitational radiation emitted by a 
cosmic string loop is still an open question. In the previ- 
ous discussion we have introduced two parameters, q, the 
spectral index and, n», the cut-off in the radiation spec- 
trum, to model the possible effects. In this section we 
investigate how these two parameters affect the observed 
spectrum. Previous works have used a range of values 
for q and which have sometimes led them in impos- 
ing very strong constraints on the string tension. For 
example, DH07 typically used n* = 1, that is, they only 
considered emission in fundamental mode of the string 
loops. Although they investigated cases up to n* s» 5, 
they found that this had only a small effect on the power 
spectrum of GWs and does not significantly effect the 
bounds on the string tension from PTAs. In contrast, 
CA92 used n* — > oo, which was done by replacing the 
summation in Eq. (|24[) with an integral in order to make 
the calculation tractable. As we will see, rx* = 1 and 
n* = oo give very different results. Damour and Vilenkin 
[5fjl \El\ and Siemens et al. [63, |62| followed a similar 
approach to CA92 making the strong assumption that 
q — 4/3 and n* = oo based on their study of cusp emis- 
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TABLE I: The percentage increase in the peak frequency and decrease in the peak amplitude of Q gw h 2 for a particular value of 
n* relative to n, = 1 and whether we are in the large (a — 10~ 4 ) or small loop (a = 10 -9 ) regimes. The results are categorized 
according to the value of the spectral index. 
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FIG. 7: Plots of fi gw /i 2 for a cosmic string network with 
Gfi/c 2 = 10~ 7 , a = 0.1 (top), a = 10~ 9 (bottom), and 
q = 4/3 for n* = 1 (black), n, = 10 4 (red), n,=5x 10 4 
(blue), n* = 10 5 (green) and n* = 10 6 (orange). In the upper 
right side of each plot a magnification of the area of interest 
is presented. 

sion. A more conservative approach was taken in [48[ 
who showed that constraints from PTAs are very sensi- 
tive to the choices of q and n» and in particular that the 
predicted spectra for q = 2 and = oo are very similar 
to q = 4/3 and n* — 10 3 — 10 4 in the nHz frequency 
range. 

As we will see, the values of ri* and q have a signif- 
icant effect on the amplitude and the slope of the GW 
spectrum in the region of radiation to matter era transi- 
tion, between the peak and the flat part of the spectrum. 



This is of critical importance for PTAs, since the fre- 
quency windows probed for the majority of Gfi/c 2 — a 
combinations of interest are in this region. The PTA fre- 
quency window lies outside this region in three cases: (i) 
Gfi/c 2 > 10~ 7 and a > 10~ 6 , where it probes the flat 
part of the spectrum, (ii) in the case of very small ten- 
sion networks, Gfi/c 2 < 10~ u independent of a, where 
it probes the region to the left of the peak (see, i.e., 
Fig. HIE]), and (hi) in the case a < 10~ 12 where the GW 
spectrum is always at higher frequencies irrespectively of 
the string tension. 

In Fig. [F5] we present the GW spectra for two repre- 
sentative scenarios varying q and large loop case, 
a = 0.1 and a small loop case, a = 1CP 9 , both with 
G/k/c 2 = 10~ 7 . The first thing to note, as found in [48| . 
is that the spectrum is relatively independent of n* when 
q = 2 both in the case of a = 0.1 and a = 10~ 9 . The 
modifications to the spectrum seen there are similar for 
the two values of a, but slightly more pronounced for 
a = 10~ 9 . Increasing n* from 1 appears to move the 
peak in the spectrum to slightly higher frequencies and 
its amplitude is also slightly reduced. For n > n sat ~ 100 
the changes in the spectrum are almost negligible, where 
we define n sa t to be the saturation point for n*, above 
which large increases of the value of n* result in negligible 
effects on the GW spectrum. 

The case of q = 4/3 is somewhat different. The shape 
of the spectra in the region of the peak is significantly 
affected by varying n*. In the case of a — 0.1 there is a 
smooth broadening of the spectrum due to loops formed 
in the matter era, with the actual peak position moving 
to higher frequencies and the amplitude being reduced. 
The situation is similar for a = 10~ 9 ; however, the spec- 
trum appears to generate a hump as n* increases. There 
does not appear to be a convergence of the spectrum 
for the values of n* presented in Fig. [5] We have inves- 
tigated the convergence of the spectrum in the case of 
large and small loops in Fig. [7J where we present spectra 
for Gfi/c 2 = 10~ 7 , a = 0.1, 10" 9 and q = 4/3 with as 
high as 10 6 . We see that in the case of large loops there 
are only minimal differences for n* > 10 4 suggesting that 
"■sat ~ 10 4 for q = 4/3. In the case of small loops a sim- 
ilar trend is observed, with the appearing hump moving 
along the radiation to matter era radiation tail of the 
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FIG. 8: The slope, dci, for cosmic string networks with n* = 10 4 (left) and n» = 1 (right) as a function of the string tension 
G/i/c 2 and the birth scale a of the loops. For both plots, q = 4/3. The dashed line corresponds to the spectral index d — — 7/6 
proposed in [58l ]. We see that for a wide range of the parameter space d is very different from —7/6. 





a > TCfifc 2 


a < TGfi/c 2 


n* 


Peak fr. 


Peak amp. 


Peak fr. 


Peak amp. 


10 


32% 


-2% 


52% 


-8% 


10 2 


135% 


-11% 


57% 


-27% 


10 3 


158% 


-19% 


57% 


-36% 


10 4 


161% 


-22% 


57% 


-36% 



TABLE II: The percentage increase in the peak frequency and 
decrease in the peak amplitude of ft gw h 2 for networks with 
q = 4/3 relative to the values for q = 2 as a function of n*. 



spectrum towards higher frequencies until it reaches the 
fiat part of the spectrum where it disappears (something 
that can be seen in Fig. [5]) . Whereas in general the spec- 
trum has converged for n* > 10 4 , the area near the tail 
of the spectrum still evolves until > 10 6 . 

In order to re-enforce the results of Fig. El in Table|T]we 
present the percentage differences in the peak frequency 
and amplitude between models with the same values of 
G/i/c 2 = 1(T 7 , a = 0.1, 1CT 9 and q = 4/3 but differ- 
ent values of n* relative to that for n* = 1. In Table ITT1 
we present similar information for a change in q, keep- 
ing all parameters fixed and changing q from 2 to 4/3. 
The results vary numerically for other combinations of 
G/i/c 2 — a, but exhibit exactly the same trend in both 
the large and the small loop regimes. 

When we discuss the constraints on G/i/c 2 due to 
observations we will be interested in the amplitude of 



the spectrum at the appropriate frequency and also the 
slope of the spectrum. If we define d and d,Q such that 
hgw(f) oc f d and fi gw cx f d ° then d n = 2(d+ 1). It is 
often suggested [58| that d — —7/6 for cosmic strings, 
but from a cursory examination of the spectra presented 
in Fig. HI El EJit is clear that this is not the case if G/i/c 2 , 
a, q and are allowed to vary. We measured e?n in 
the region between 31 nHz and 32 nHz as a function of 
G/i/c 2 and a for q = 4/3 and two values for n*, n* = 1 
and n* = 10 4 . In general, realistic PTA experiments are 
sensitive to GWs with frequencies of a few nHz. How- 
ever, later in this work we will calculate the string tension 
constraints based on limits placed for a frequency (lyr) -1 
and therefore we have chosen this range. 



In Fig. [8] we present a plot of dii(Gfi/c 2 ,a) at a fre- 
quency / = (lyr) -1 for cosmic string networks with 
ri„ = 1, 10 4 and q = 4/3. The results for other values 
of n* and q are of course numerically different, but they 
follow a very similar trend. We see that describing the 
cosmic string GW spectrum with a simple power law in 
the frequencies probed by PTAs is far from reality. In the 
tension range G/i/c 2 > 10 -11 the slope is generally neg- 
ative, and PTAs probe the whole area of the spectrum 
which lies between the matter era peak and the radia- 
tion era fiat part of the spectrum. However, when we 
go to small tensions, G/i/c 2 < 10~ n , the PTA frequency 
window falls to the left of the matter era peak and the 
spectrum slope becomes positive. 
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4- Varying p 

In order to calculate the GW spectra of cosmic string 
networks with p < 1 we used the same one-scale model 
based code, only changing the value of A as in Eq. (1231) . 
In Fig. [9] we present the results for different values of 
the intercommutation probability (p = 0.1,10~ 2 ,10~ 3 ) 
and for both scaling laws (k = 1 and k = 0.6). In all 
computations we used G/x/c 2 = 10~ 7 , a = 0.1, n.* = 1 
and q = 4/3. 

The reduced intercommutation probability leads to an 
increased number density of cosmic string loops, and 
therefore, an increased number of GW sources which give 
higher SGWB. The uniform scaling across the frequency 
band make the effects of decreasing p similar to those of 
increasing G/i/c 2 in the small loop regime, see Fig. [4] 



a to vary as an unknown parameter. 

There is, of course, nothing to prevent loops being born 
over a ran ge of diffe rent scales, both small and large. Re- 
cent work jl03lll04l | has presented arguments supporting 
the idea that loops are created at two different scales, 
with 90% of the loops created at the gravitational back- 
reaction scale (a ~ TG/z/c 2 ) and 10% at large scales 
(a « Q.l). In this section we will discuss the effects of 
relaxing the assumption of a single loop production scale. 
We will consider two possibilities: the firs t is a two -scale 
model for loop production motivated by [103l Il04j , and 
the second is for the initial loop distribution to have a 
log- normal distribution with the mean being oc t^. 



1. 2-scale networks 




log 10 f (Hz) 

FIG. 9: The effects of varying the intercommutation proba- 
bility p for Gfi/c 2 — 10 -7 and a = 0.1. The red, green and 
blue lines are for p — 0.1, 10 -2 , 10 -3 respectively for k — 0.6. 
With the same color scheme but with dashed lines we show 
the equivalent results for k = 1. The black solid line is for 
p=l. 



C. Improved modeling of loop production 



The one-scale model can be easily adapted in order to 
describe networks with more than one scale for the new- 
born loops. A network that produces loops of two scales, 
will create the same shape of SGWB as the one created 
by two distinct one-scale networks, but with the relative 
amplitude of the two carefully normalized. In order to 
impose this normalization we need to enforce the condi- 
tion that a specific amount of energy has to be "lost" 
from the network in the form of loops per unit time; a 
requirement necessary for scaling. If we want to describe 
the loop distribution with two separate networks of a.\ 
and ai respectively, the amount of energy transferred to 
loops per unit time will be 



dE- 



2— scale 



dt 



dEx 



,dE- 



^i + c^, (31) 



dt 



dt 



where dEi/dt, dE^jdi are the energies transferred to the 
two different types of loops and cf , cf are the appropriate 
weighting factors. The values of cf,C2 are the relative 
amounts of energy we want to be channeled in the loops 
of each scale. 

From Eq. (0 we can see that the energy channeled in 
loops of a specific scale is proportional to the number of 
loops created 



It is extremely difficult to model the distribution of 
loops produced by a cosmic string network. The cal- 
culations presented so far are based on the assumption 
that loops are produced with a single size relative to 
the horizon, a. In the initial work on cosmic strings 
it was believed that large loops are born, with a 0.1 
[661 I90T ] . The work that followed [74| [9l[ argued against 
this, supporting the idea of a smaller length scale close to 
the scale of the gravitational backreaction a « TG/i/c 2 . 
Subsequent work has led to a somewhat confusing sit- 
uation. Some appear to sugges t lar ge loops |9a - l9a | 
(a fa lO^-lO -3 ) while others [97l - ll0l| support the view 
that the loops are small (a < TGfi/c 2 ). S ome even sug- 
gest microscopic loops with ss 5 [ll], [jj Il02j | , where 5 
is the string core width. For this reason we have allowed 



cLEi.2 diVi,2 

oc 



dt 



dt 



(32) 



and therefore the normalizing constants cf 2 can a ls° be 
used as normalizing constants for the loop number den- 
sity 



'1,2 



'1,2 - 



(33) 



In this way we can calculate the SGWB for each network 
individually, using the method described before, but this 
time we will normalize the number density of loops with 
and c 2 respectively. Once we have calculated the 
individual SGWBs, we add them and we get the SGWB 
produced by the 2-scale network. 
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log 10 f (Hz) log 10 f (Hz) 

FIG. 10: Qgwh 2 for networks which produce two different sizes of loops corresponding to scales ot\ and a?2 relative to the horizon. 
In all cases we used G/i/c 2 = 10~ 7 , q = 4/3, n» = 1 and p = 1. The ax, Q.2 and the relative percentages used are shown in 
the individual figures. The black lines are the total spectra for the overall 2-scale network. The solid blue and red lines are 
the individual contributions to the total spectrum from loops with initial sizes relative to the horizon, Qi and «2 respectively. 
For reference we have also included spectra for the networks with ot\ and a?2 (blue and red dashed lines respectively) assuming 
that all of the energy was channeled into each of the individual loop sizes. 



In Fig. \W\ we present results for representative mod- 
els for which newborn loops are created with two dis- 
tinct different scales, a\ and GV2 and all other parame- 
ters given by their fiducial values. The solid blue and 
red lines are the individual contributions due to loops 
with ot\ and ai respectively, normalized to their relative 
contributions to the total emission. For reference, the 
dashed lines are the results we would expect if all the 
energy was channeled into loops with either a\ or 0.2- 
In Fig. \TU[ i) the spectrum corresponds to a case with 
ai,Q!2 > TGfi/c 2 . In this case, even for a very small 
amount of energy (10%) channeled into the creation of 
the largest scale loops (ai = 0.1), the GWs emitted by 
them dominate the overall result, with the smaller scale 
loops (a2 = 10~ 4 , 90%) only dominating the overall spec- 
trum at very low frequencies, / < 10~ 12 Hz. This effect is 
even more obvious in Fig. HOf iii). where there is a larger 
percentage of the energy channeled into the creation of 
the large scale loops (a>i — 10~ 2 , 40%). There we see that 
the contribution from the large loops contributes most of 
the overall result. From these results, and the other cases 
we have investigated but not included in the figure, we 
conclude that the spectrum of the 2-scale model has a 
similar shape to ones created in cases when there is only 



one loop production size when ai, a-i > TGfi/c 2 . 

The situation is very different if a\ > TGfi/c 2 > CV2, 
as illustrated in Fig. [TOT ii). In that case, we find that 
the spectrum of the SGWB has two peaks, one due to 
each loop population. This is consistent with the discus- 
sion of Sec. IIII B 21 where we showed that the peak of 
the spectrum for a < TGfi/c 2 is shifted towards higher 
frequencies. This means that its peak will be higher than 
that due to the emission from the loops with a > TGfi/c 2 
in a frequency higher than the / pca k of the a > TG^l/c 2 
loops. In Fig. fTOlf ii). 90% of the energy goes to the cre- 
ation of small scale loops (02) and therefore the peak 
corresponding to the a-2 (higher frequency peak) will be 
more prominent. In general, which of the two peaks dom- 
inates over the other depends on the amount of energy 
that the network allocates to the respective loop sizes. 
The same behavior is exhibited in Fig. [TUf iv). where 
cti,a>2 < TGfi/c 2 . We see the two peaks are even more 
clear in this case, with the second being slightly higher 
since a larger amount of energy goes into the creation of 
loops at that scale (02)- Based on a more general inves- 
tigation of the parameter space, the basic feature of two 
peaks is something which appears to be generic in the 
case when one, or both, of the loop production scales is 
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Gfi/c 2 = 1(T 7 , a, = 0.1. a 2 = 5x10"' 
G/i/c 2 = 10" 8 , a, = 0.1, a 2 = 5 x 10~ 7 
/c 2 = 10"', ff, = 0.1, «2 = 5 x 10~ 8 
G/i/c 2 = 10~' , a, =0.1, 0-2 = 5x10-' 



then the overall spectrum is 



log 10 f (Hz) 



FIG. 11: Qg W h plots for 2-scale networks with 10% of the 
energy lost into the production of large loops (qi = 0.1) and 
90% to loops at the gravitational backreaction scale (a?2 = 
TGfj,/c 2 ). GW spectra for various Gfj,/c 2 values are presented. 
As we go to smaller tensions, the sharp peak created by the 
small loops starts to weaken until we reach the value G/1/c 2 — 
10 -10 , below which the total GW spectrum of the 2-scale 
network is actually created only by the large loops. 



lower than the gravitational backreaction scale. 

In Fig.[TT]we present results for networ ks which specif- 
ically follow the predictions of [l03l Il04| . in which 10% 
of the energy lost by the network goes into large loops 
of scale ai « 0.1 and the remaining 90% goes into loops 
at the gravitational backreaction scale (a2 « TG/x/c 2 ). 
We present results for various values of G/i/c 2 . For 
G/x/c 2 > 10~ 10 , we find a two-peak spectrum, as in 
Fig. HOt ii). but much less prominent. The position of the 
peak created by the small loops is almost the same as the 
one created by the large loops, creating a spectrum simi- 
lar to that of a pure one-scale network but with a sharper 
peak. However, the peak amplitude continuously de- 
creases as we decrease G/z/c 2 , and when Gfi/c 2 < 10~ 10 
the peak disappears and the large scale loops' contribu- 
tion totally dominates the spectrum. 



2. Log-normal distribution for loop production 

Another possibility to consider is that the loops are 
born with size £(t^) = f T ad^(t\,) but with a having a 
distribution, 'P(cv), an idea qualitat ively justified by some 
recent simulations [75l IfM loBL Hfj3| . 

One can model any distribution using an adaptation of 
the methods used in the earlier sections by splitting the 
distribution into N populations with loop production size 
relative to the horizon at for i = 1, .., N and a fraction of 
loops in each bin of size Aa given by Vipii) Aa, assuming 
that 



poo 

/ daV(a) = l. 
Jo 



(34) 



nsw(/) 



N 



(35) 



where the c$ = V(a,i)Aa are computed in order to enforce 
the overall energy loss required to maintain scaling. 

We have chosen to use the log-normal distribution to 
model the loop distribution 



V(a) = 



1 



g v 2ttoi 



exp 



1 


log 




2" 




(«o)_ 




~2a* 







(36) 



If f2gw (/) is the spectrum computed from loops with a 



which corresponds to a Gaussian distribution in log a 
with mean log ao and variance a. Typically, we will dis- 
cretise the distribution in the range —9 < log 10 a < —1) 
with 21 bins with A(log 10 a) = 8/21. We also exper- 
imented with smaller bin sizes, but we didn't see any 
change in the final results, making our choice the most 
computationally effective. 

The results we present here use the fiducial set of pa- 
rameters, except that we consider the two cases n* = 1 
and n* = 10 5 , and we have varied cvo and o. To ex- 
hibit the important effects, we chose to use three differ- 
ent mean values, ao — 10~ 3 (a > TGfi/c 2 ), a = 10 -8 
(a < TGfi/c 2 ) and a = 5 x 10~ 6 {a = TG/i/c 2 ), and 
two values of the variance a = 0.4 and a = 0.8 to demon- 
strate the behavior of multiple scale networks in the case 
of large loops, small loops and those produced at the 
gravitational backreaction scale respectively. 

In Fig.[T2]we present the results of our computations. 
In the upper right of each plot is the mean value of a 
used. The red lines are for the power spectra of a one- 
scale model whose value of a equal to the mean used for 
the multiscale networks. The green and blue lines are 
the spectra for the multiscale networks with a = 0.4 and 
a = 0.8 respectively. The spectra have different behavior 
in the case of large and small loops. In the case of large 
loop creation (ao > TG/x/c 2 ), the multiscale networks 
give higher SGWB amplitudes than the corresponding 
one-scale networks and moreover, the increase is higher 
for larger values of a since the spectrum is dominated 
by the few very large loops. This amplitude increase is 
seen over the whole frequency region for / > / pca k- This 
behavior reverses in the case of networks with small loops 
or loops near the gravitational backreaction scale (ao < 
TGfi/c 2 ). The flat part of the spectrum, as expected, 
is not affected when all the loops are born small (a < 
TG/i/c 2 ). However, when «o = TGix/c 2 a small increase 
in the amplitude of the flat part is seen due to the loops 
at the tail of the distribution, born with a > TG/i/c 2 . 
Although there are very few of them, they can have a 
noticeable effect on the spectrum. 

On the other hand, the peak region created from GW 
emission during the matter era exhibits a richer behavior. 
A common feature is the decrease of the peak amplitude 
combined with a broadening of the whole peak region, 
when compared to the corresponding results from the 
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FIG. 12: f2 gw /i 2 plots for networks in which cosmic string loops are born with a given by a log-normal distribution. All the 
computations were performed for a network with G/x/c 2 — 10~ 7 , q = 4/3, n* = 1 and p = 1. In the left column are the plots 
with n* = 1 and in the right column are results with n„ = 10 s . The green and blue lines are created using the log- normal loop 
production distribution with a mean value of ao and with a = 0.4 and a — 0.8 respectively. The red lines are the results for a 
single loop production size a — ao. 



one-scale model. This is even more prominent in the 
case where ao < TG^/c 2 , because of the shift in the 
peak frequency from the loops born with different values 
of a. Interestingly, we see that the low frequency cut-off 
of the one-scale network is no longer present, with the 
SGWB spectrum extending to lower frequencies. This 
extra emission into these low frequencies is created by 
the few large loops at the tail of the distribution which 
individually, have lower cut-offs. Finally, as seen from 
Fig. [T21 the number of emission modes does not have 
any significant effect on the spectra of the multiscale net- 



works but in general suppresses the effects of amplitude 
decrease of the flat portion of the spectrum, peak am- 
plitude decrease and peak region broadening already dis- 
cussed. 



IV. PULSAR TIMING CONSTRAINTS ON THE 
COSMIC STRING TENSION 

In this section we will use the constraints on the SGWB 
from pulsars to impose a constraint on the dimensionless 



Authors 


Q gvi h 2 bound 


Bertotti, Carr, Rees [107] 


1.0 x 10~ 3 @20nHz 


Kaspi, Taylor, Ryba [108] 


6.0 x 10 _8 @4.5nHz 


Thorsett, Dewey [106] 


1.0 x 10~ 8 @4.5nHz 


McHugh, Zalamansky, Vertotte [109] 


9.3 x 10" 8 @4.5nHz 


Lommen [110] 


2.0 x 10 _9 @1.9nHz 


Jenet et al. [61] (PPTA) 


1.9 x 10" 8 @4.0nHz 


van Haasteren et al. [65] 


5.6 x 10~ 9 @4.0nHz 



TABLE III: The f2 gw (/)ft 2 limits published in the literature 
so far. The Jenet et al. limit is derived for supermassive black 
holes (dn = 2/3) and cosmic strings (dn = —1/3) as sources of 
the SGWB. The van Haasteren et al. limit quoted here, is for 
supermassive black holes at a frequency of 4.0nHz, for direct 
comparison with the Jenet et al. limit. The Thorsett and 
Dewey limit was criticized by McHugh et al. because of the 
statistical approach used. The Lommen limit, although the 
most stringent so far, was also criticized for similar reasons. 

cosmic string tension parameter G/i/c 2 using the EPTA 
SGWB limit [65[ and discuss the p ossibi lities of using 
the future constraints from the LEAP [l05j j project which 
should be available in the near future. There is a lon g lit- 
erature on this subject [H, M, 13 M, MM, M, M,UM 
based on many different assumptions for the string net- 
work evolution, radiation emission, and pulsar timing 
data. Often there are contradictory constraints based on 
the same data since some authors make very strong as- 
sumptions about the expected SGWB spectrum whereas 
others are more conservative. This is at best confusing to 
those uninitiated in the details of string evolution, some- 
thing which this paper is aimed at clarifying. Therefore, 
we will take the view that none of the parameters (a, q 
and n») which we have described in the earlier sections 
are known, extending the work done in |48| . The head- 
line constraints which we will quote will be the highest 
upper bound possible for any parameter, this being the 
absolute constraint on G\ij(? - one that is conservative 
and which we can never go back on! In addition we will 
present constraints for various specific models which have 
been discussed in the literature. 

A. EPTA constraint on O gw /i 2 

The SGWB limit we will use is that which comes from 
the EPTA who have computed a 95% exclusion curve 
for the SGWB amplitude h gw as a function of the local 
slope, d, at / = (lyr) -1 « 32nHz. Previous work on the 
subject has typically quoted an upper limit on Q gw h 2 at 
a particular frequency. Over the last 30 years various 
data from PTA experiments have set constraints, as pre- 
sented in Table IIIII The limits which we will discuss in 
the subsequent sections will use the full exclusion curve 
from [65f . However, we shall follow a different approach 
in the limit predictions for the LEAP project since such 
an exclusion curve is obviously not available yet. In order 



17 

to have a conservative projection for LEAP, we will use 
a spectral index dn = to make a projected bound. By 
making this choice we guarantee that the constraints on 
n gw (/) will be the most conservative applicable to cos- 
mic strings for the majority of G/i/c 2 — a combinations. 
Note that dn < in the radiation to matter transition 
epoch frequency range for any cosmic string model. We 
will come again to the applicability of this idea when 
we will calculate the projected LEAP constraints. More- 
over, the LEAP constraints are evaluated at a frequency 
of (5yr) _1 , which is calculated from the duration of a 
typical PTA experiment. 

In [65| the strain of the SGWB is described by a power 
law of the form h gw — A(f /yr~ 1 ) d and an upper bound 
A = 6 x 10~ 15 was established in the case of supermassive 
black holes (d = —2/3) at / = (lyr) -1 . It is interest- 
ing to estimate the corresponding constraint on fl gw h 2 
in order to confirm that the EPTA SGWB limit is the 
strongest possible one with which to per form our analy- 
sis. The previously published limit in [61| (see Table IlLi]) . 
was quoted for d = —2/3 at / = (8yr) _1 . For the same 
parameter values, the projection of the EPTA limit gives 

tt g w,EPTA(/)/i 2 £ 5-6 x 1CT 9 , (37) 

at 95% confidence level, a significant improvement in 
comparison to the Jenet et al. limit [61]. 

We note that there is a constraint on G^t/c 2 presented 
in 

G/Vc 2 < 4.0 x 10" 9 , (38) 

that claims to improve the previous one by Jenet et al. 
[6ll ] which was 

Gfi/c 2 < 1.5 x 10~ 8 . 

However, both of these limits arc based on the approach 
of Damour and Vilenkin [58| who make strong assump- 
tions about the amplitude and slope of the SGWB from 
strings. In particular, they describe the amplitude of the 
cosmic string loop generated SGWB using a model based 
on the emission from cusps on cosmic string loops with a 
slope O gw /i 2 oc /~ 7 / 6 . As we have already discussed, loop 
decay solely through cusp emission is a rather strong as- 
sumption, and the possible behavior of the cosmic string 
GW spectrum is much richer than this. In addition, the 
analytic results of [58[ are based on a simplified cos- 
mic string evolution model (i.e., the loop number den- 
sity n(£,t) calculation). Damour and Vilenkin assume 
n(t) ~ (TGiit 3 /c 2 )" 1 independent of a, whereas the an- 
alytic model of [63, [68| which is compatible with our 
computations yields n(t) ~ (1 + TGfia^ 1 c^ 2 )(at 3 )^ 1 . 
This agrees with that of Damour and Vilenkin when 
a = TG/i/c 2 , that is, the loop production is at the grav- 
itational backreaction scale. 
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value at a frequency / = 
= 4/3 and n* = 1 (black), 
= 10 2 (orange). The purple dashed curve is the analytic approximation 



FIG. 13: 

(lyr) -1 in the G/i/c 2 — a parameter space. The solid lines are for the EPTA (lyr) -1 limit with q 
n* = 10 3 (red), n, = 10 4 (green) and for q = 2, 

of the fiat part of the spectrum as described in [l2j, again for the EPTA limit on S7 gw /i 2 . The dot-dashed brown line is the 
present LIGO limit. The constraint set by LIGO is almost independent of n* and q and is also only very weakly dependent on 
a. The dashed lines are the corresponding curves the planned LEAP sensitivity at a frequency / = (5yr) -1 . The thick dashed 
black line shows a — FG/i/c 2 . 



B. Conservative constraint on G/i/c 2 

In order to calculate the constraints on the cosmic 
string tension, we have to find all the cosmic string net- 
work parameter combinations which lead to a SGWB 
amplitude which is in agreement with the EPTA SGWB 
limit at a frequency of (lyr) -1 . For this, we computed 
the SGWB (Q ew h 2 (f)) for more than 3000 parameter 
sets covering almost all the range of the theoretically ex- 
pected values for G/i/c 2 , a, n*,q and p. For each set 
of fixed n* , q and p, we deduced from our computations 
the quantity f2 gw (G/z/c 2 , a, dn)h 2 , which gives the am- 
plitude of the SGWB as a function of Gfi/c 2 ,a and the 
local slope da at a frequency of (1 yr) -1 . The EPTA limit 



is given in the form of a 95% exclusion curve of the form 
^gw.EPTA(rf)- From this, we calculated the correspond- 
ing f2 gw ,EPTA(rffi)^ 2 95% exclusion curve which would 
be applied to our results. The G/z/c 2 — a combinations 
which provide the constraint curve for each set of n*, q 
and p in the G/i/c — a parameter space was calculated 
by requiring fi gw (G^/c 2 , a, d n )h 2 = £l gw ,EPTA(dn)h 2 . 

In Fig. [T3] we set p = 1 and present constraints for 
various values of q and which satisfy the EPTA limit. 
Specifically, we have plotted the cosmic string network 
families for n» = 1 (blue); n* = 10 2 and q = 4/3 (red), 
n» = 10 3 and q = 4/3 (green); n* = 10 4 and q = 4/3 
(black); rt* = 10 2 and q = 2 (orange). The dashed 
black line is a = TGfi/c 2 separating the large and the 
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small loop production regions. The constraints for high- 
n* models are stronger than those with low-n* for most 
of the a — Gp/c 2 combinations except for the case of very 
small G/i/c 2 or a, where the opposite takes place. In the 
same figure we also present (purple dashed line) the con- 
straints of the analytic approximation presented in (l2l | 
that is just from the radiation era contribution. This is 
a good approximation to our results in the case of large 
loops, however it is only a conservative upper limit for 
lower values of a. 

The most conservative and generic constraint on the 
cosmic string tension can be set by the curves presented 
in Fig. 1131 This is provided by the cosmic string networks 
with a ~ 10~ 5 and n* = 1, and is 

Gp/c 2 < 5.3 x 1CT 7 , (39) 

which is a 95% upper bound for this specific set of pa- 
rameters. 

For comparison, we have also used the bound on the 
SGWB set by LIGO to obtain a constraint on the string 
tension. We expect this to be worse than the one set 
by PTAs not only because the limit on fi gw /i 2 itself 
is smaller, but also because the LIGO frequency band 
(~ 1kHz) is probing the radiation era part of the spec- 
trum which has a smaller ampl itude than the matter era 
equivalent. The LIGO limit is [HO] 

^gwxico^ 2 < 3.6 x 1CT 6 , (40) 

at 95% confidence level. In Fig. [13] with the dot-dashed 
brown line we show all the cosmic string configurations 
compatible with this bound at 1kHz. For p = 1 the most 
conservative constraint is 

G/i/c 2 < 2.6 x 1CT 4 . (41) 

Advanced LIGO is expected to make sufficient improve- 
ment, although it is unlikely to compete with PTAs in 
the short term. Conversely, one of its biggest advantages 
is that at these frequencies the spectrum is independent 
of the modeling of the string radiation spectrum, that is, 
q and n*. 

The prospects of improving the PTA limit in the near 
future are very promising. As part of the EPTA, the 
LEAP [l05j | project, is a large collaboration of the major 
telescopes of the EPTA members (Effelsberg Telescope, 
Lovell Telescope, Nancay Telescope, Westerbork Synthe- 
sis Radio Telescope, Sardinia Radio Telescope) which will 
use long-baseline array techniques to provide pulsar tim- 
ing data equivalent to that of a 200-meter telescope. It 
is estimated that LEAP data will push t he s ensitivity 
in GWs down to fl gw , leap {f)h 2 < 10" 10 OH in a few 
years. In Fig. \W[ with dashed lines and the same color 
scheme as for the EPTA limits, we present the constraints 
one might expect from such a limit on the SGWB. If we 
assume that LEAP does not make a detection, then we 
would expect to set a limit of G/i/c 2 < 2.0 x 10~ 9 ; an 
improvement of more than two orders of magnitude. The 
careful reader will notice that for the LEAP constraint 



the order of the lines is reversed for high-n* values and 
the high-n* models are above the low-ro* cases. This is 
expected, since we are in the low a — Gp/c region of the 
parameter space and the PTA frequency window probes 
the part of the spectrum which is on the left side of the 
peak, where we can see from Fig. [5] that configurations 
with high-rt* give lower amplitudes. 

We will conclude this section with some comments on 
our choice of a spectral index e?n = for the calculation 
of the projected LEAP constraints. As we mentioned 
in Sec. IIII B 31 the PTA window probes the area on the 
left of the peak (positive slope) for cosmic string net- 
works with G/i/c 2 < 10~ u . The constraint lines from 
the projected LEAP limit correspond to models with 
G/i/c 2 > 10~ u in the range a € [10~ 9 ,10~ 3 ]. For all 
these models, indeed our assumption behaves well and 
moreover, it gives slightly overestimated values for G/i/c 2 
since the actual slope in this region is negative. The pro- 
jected tension constraint is given by the cosmic string 
networks with n* — 1 (blue dashed line in Fig. I13[) which 
lie in this region, so our results are robust. On the other 
hand, in the region a € [10~ 3 , 0.1] the cosmic string mod- 
els have G/t/c 2 < 10~ n and their slope is positive. For 
these networks, our assumption is invalid and the ten- 
sions are slightly underestimated. This under-estimation, 
however, is not sufficient to overcome the much stronger 
constraints provided by the n* = 1 networks. 

C. Conservative tension constraints for p / 1 

We have used the EPTA h gw 95% exclusion curve to 
calculate conservative constraints on G/i/c 2 for the cases 
of p < 1 for both scaling laws, poo oc P~ l and p^ oc p~ 6 
and the results are presented in Fig. fT?] In this plot, 
the curves correspond to the highest tension values only, 
meaning that a significant part of each curve is provided 
by the n* = 1 case (in approximately the mid-a range) 
and the rest by = 10 4 (in the low and high a ranges). 
The constraints for the p^ = p~ Q - 6 scaling law are plot- 
ted with solid red, blue and green lines for the cases of 
p = 0.1, 10 ~ 2 and 10~ 3 respectively. The dashed lines 
are the corresponding results for the p^ = p^ 1 scaling 
law. The conservative G/i/c 2 constraints are presented in 
Table IIVI The results are indicative of how sensitive the 
tension constraints are, not only for the value of the in- 
tercommutation probability, but also to the exact scaling 
law as well. 



p 


G/i/c 2 


G/i/c 2 




( Poo =p-°-») 


(Poo =p~ v ) 


0.1 


1.2 x 10~ 7 


4.4 x 10" 8 


0.01 


2.0 x 10~ 8 


1.3 x 10" 9 


0.001 


2.8 x 10" 9 


9.3 x 10~ 12 



TABLE IV: Cosmic string tension constraints for models with 
p < 1 and for the two scaling laws discussed in the literature. 
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FIG. 14: The exclusion limits for cosmic string networks with different intercommutation probabilities and different scaling 
predictions (p^ oc p _1 and p^ oc p~ ' 6 ). These curves where created by combining the equivalent lines for n* = 1 (for most of 
the range of a) and n* = 10 4 (for large/small a). All of these network configurations give an amplitude equal to S! gw ,EPTA/i 2 
at the frequency of (lyr) -1 . The solid black line is for networks with p = 1 and the solid red, blue and green are for networks 
with p = 0.1, 10 -2 , 10 -3 respectively for poo oc p~ 6 . The dashed lines are the networks with p x oc using the same color 
scheme. The thick dashed black line shows a — YGp,/c 2 . 



D. Constraints on specific scenarios 

In addition to the previous, generic constraint on the 
string tension, we have also computed constraints for 
some specific cosmic string network scenarios popular in 
the literature. These models are based either on evo- 
lution specific simulations of cosmic string networks or 
theoretical arguments which estimate the birth scale of 
the loops, a. Unfortunately, as we mentioned before no 
definite conclusion has been drawn yet. Depending on 
the value of a we can separate these models into four 
categories 

i. Large loops: it has been suggested that most of the 
cosmic string loops are born with a size compara- 



ble to the horizon radius, with a « 0.1 [95ll96{. Such 
simulations have been performed for both flat and ex- 
panding spacetimes. The production of small loops 
along with the large ones has also been observed in 
these simulations but it was suggested that they are a 
transient effect which disappears when the dynamic 
range of the simulations is increased. The most up- 
to-date simulations in this category are those in |75| 
where a slightly different value of a s» 0.05 was sug- 
gested. 

ii. Intermediate loops: in some simulations [92], |93[ it 
has been argued against the creation of such large 
scale loops supporting the idea that the smaller loops 
are not a transient phenomenon but instead domi- 
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nate the loop production. Both of these simulations 
were performed for expanding spacetimes and they 
concluded that a ps 10 ~ 2 — 1CP 3 without observing 
any significant production of large scale loops. 



iii. Loop production is governed by the gravitational 
backreaction scale : in this scenario cosmic string 
loops are created at the scale of gravitational back- 
reaction a ps TGfj,/c 2 [74| as a direct result of the 
small-scale structure on the cosmic string network. 
In more recent simulations |98j it has been proposed 
that the scale is even smaller and isax (TGjj,/c 2 ) k 
with k = 3/2 in the radiation era and k = 5/2 in the 
matter era. 



Loop scale (a) 


Gfi/c bound 


Ref. 


0.1 


6.5 x 10" 11 


[95, 96] 


0.05 


8.8 x 10 -11 


[75] 


10~ 2 


7.0 x HT 10 


[92, 93] 


TGfi/<? 


5.3 x 10" 7 


[74] 


10% o; — 0.1 + 90% a = TGfi/c 2 


4.1 x 10" 8 


[103, 10J 


~10" 9 


1.9 x 10" 8 





TABLE V: Cosmic string tension limits for individual cosmic 
string models with specific a values predicted by particular 
simulations. 



V. CONCLUSIONS 



iv. Two scale loops: Some simulations [l03l Il04j sug- 
gest that cosmic string loops are born with two dif- 
ferent scales. Specifically, 10% of the string energy 
converted to loops per unit time is channeled into 
loops with a ps 0.1 and the rest 90% into loops with 
a ps TG/i/c 2 scale loops. 



There are also simulations which suggest that cosmic 
string loops are born at a fixed size that does not scale 
with cosmic time and it is equal to the string width, that 
is ol ps 0. Such tiny loops are expected to decay mainly 
through particle emission [72j and are of no interest for 
this work. However, even if their GW emission is signif- 
icant, it will take place at frequencies much higher than 
those probed by PTAs. 

We have computed the GW spectra of cosmic string 
models in each of these scenarios in the tension range 
Gii/c 2 G [10" 5 ,10- 16 ] and applied the EPTA 2<r-limit 
on the SGWB at the frequency of /epta = (lyr) -1 to 
set constraints on the string tension. The results for all 
categories are presented in Table IVl Constraints for mod- 
els with a ps (TGn/c 2 ) k have not been computed since 
their low frequency cut-off is at frequencies higher than 
those probed by PTAs, rendering them unobservable. In 
the categories (ii), (iii) and (iv), the most conservative 
constraint is given by n* = 1, whereas in case (i) it is 
given by = 10 4 (sec in Fig. Q21 that these networks 
give higher tensions near a = 0.1). 

In the case of large loop production which is favored 
by the most recent simulations, a very stringent limit of 
G/z/c 2 < 6.5 x 10~ n is obtained. This limit is slightly 
weaker than the corresponding one found by the authors 
in 12], G/i/c 2 < 5.0 x 10 -11 which is based on the ana- 
lytic approximation of the amplitude of the radiation era 
part of the spectrum. In the regime a -C TGfi/c 2 , the 
smallest loops that PTAs are sensitive to have a ~ 10~ 9 . 
From Fig. fT3l we get a constraint of G/x/c 2 < 1.9 x 10~ 8 
for such loops. 



The constraints on the cosmic string tension from 
PTAs suffer from the many uncertainties concerning the 
parameters describing the loop production size, number 
density and GW emission properties. Particular treat- 
ments of the emission mechanism presented in the litera- 
ture has led to different conflictin g co nstraints [60l l63l . H3| 
with the most stringent of them [6l[ l6o| being based on 
rather strong assumptions [56l - [58l |. In this paper, we have 
expand on the work performed in (48|, performing a de- 
tailed investigation of the constraints for a wide range 
of scenarios. Using the recent limit on the SGWB at a 
frequency / = (lyr)" 1 set by the EPTA [65j, we have 
managed to set an absolute, model independent upper 
limit on the string tension without making any assump- 
tions for the emission properties. This limit is 

G[i/c 2 < 5.3 x 10~ 7 , 

at 95% confidence level. Such an approach is particularly 
necessary in the case of PTAs, since the frequency range 
they probe is potentially sensitive to all the network and 
emission parameters. Additionally, we have calculated 
constraints for cosmic string networks with p < 1 and for 
networks with specific loop birth scales which have been 
proposed in the literature. 

In order to achieve this, we investigated the effects on 
the SGWB for each of the cosmic string model parame- 
ters in the range of values which are interesting for PTAs. 
This allowed us to delineate the fundamentally different 
SGWB behavior of cosmic string networks in the large 
and small loop production regimes, defined by the grav- 
itational backreaction scale T Gfi/c 2 . Of special interest 
to the PTAs is the low frequency cut-off, which can ren- 
der the GW signatures of small loop size cosmic string 
networks undetectable even for high tension values. Ad- 
ditionally, we considered extensions to the standard one- 
scale model in order to investigate the differences in the 
SGWB created by more realistic cosmic string networks 
with a 2-scale and a log-normal loop production function. 

It is worth comparing our limits to those placed on 
G/z/c 2 from other observational approaches. 
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CMB observations: cosmic strings could be respon- 
sible for a fraction of the CMB anisotropics due 
to the Kaiser-Stebbins effect [3l[. Data from vari- 
ous CMB experiments have been used to place con- 
straints on the string tension, the most stringent of 
which arc those set by combining 7-year Wilkinson 
Microwave Anisotropy Probe (WMAP) data and 
Atac ama Cosmology Telescope (ACT) observations 

EH 

Gfi/c 2 < 1.6 x 1CT 7 , 

at 95% confidence level. This uses the methods de- 
scribed in [IJJ designed to model the simulations 
of [qTJ] . Weaker constraints are claimed in the case 
of Abelian-Higgs simulations 



114], the most con- 



straining of which is Gfi < 4.2 x 10 -7 , derived using 
all available CMB data. Sec 12] for a discussion of 
the veracity of these different approaches. 

• Gravitational lensing: cosmic strings can cre- 
ate gravitational lensing events due to the coni- 
cal shape that spacetime acquires globally around 
them P, [28[. Such events have special features 
which can distinguish them from "normal" lenses 
(that is, no magnification, identical images, odd 
number of images) and can be observed in searches 
for gravitational lenses. The m ost recent constraint 
set from such searches was in [115L 1 1 1 61 ] where the 
authors used archival data from the HST (GOODS 
and COSMOS surveys respectively) to search for 
lensing events that could have been caused by 
straight cosmic strings, managing to set a bound 
slightly weaker than the CMB one 

Gfi/c 2 < 3.0 x 1CT 7 , 

at 95% confidence level. We note, that these con- 
straints are also sensitive to some of the modeling 
issues which are important in the CMB case. 

So far, these constraints were considered more reliable 
than those from PTAs, due to the many uncertainties 
concerning the cosmic string loop distribution and the 
GW emission assumptions of previous implementations. 
In this paper we believe that we have overcome all these 
difficulties and provide equally reliable, PTA constraints. 
In addition, the prospects for PTAs are much better than 



those of other experimental approaches. In the very near 
future, LEAP is expected to improve on the present up- 
per limit on Gfi/c 2 by more than two orders of mag- 
nitude to Gfi/c 2 < 2.0 x 10~ 9 , if no detection of the 
SGWB takes place. This is a stronger constraint than 
any other expected from present or near future CM B ex- 
periments. The projecte d lim it for Planck satellite |l!7l ] 
is Gfi/c 2 < 6.5 x 10~ 8 [HU, approximately one orde r 
of magnitude better than the present one. In [Mll20l. 
the authors come to similar pessimistic conclusions about 
the prospects of present and future planned CMB exper- 
iments to detect CMB polarization. 

Future PTA experiments could provide an unprece- 
dented insight on the physics of cosmic strings. It is clear 
that projects like LEAP and the growing International 
Pulsar Timing Array [12l| will likely make a detection 
of the SGWB from supermassive black hole binaries on 
the timescale of the next five years. However, the sen- 
sitivity that will be achi eved with the next generation 
radio telescope, the SKA [l22j , will give us the opportu- 
nity to study the GW spectrum in detail. As it has been 
shown in this pap er for cosmic strings and for supermas- 
sive black holes in [l23l | , measuring the shape of the spec- 
trum is essential to be able to both distinguish between 
the source of the SGWB and also to extract more infor- 
mation about the source. In the case of supermassive 
black hole binaries, information about the supermassive 
black hole population and the evolution of binary sys- 
tems can be extracted. In the case of cosmic strings, we 
could potentially determine all the fundamental model 
parameters (loop birth scale a, intercommutation proba- 
bility p) and get a definitive answer about the exact GW 
emission mechanism and small scale structure of cosmic 
string loops (i.e. cusps or kinks). 
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